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INTRODUCTION 


!  ^ 


The  phenomenon  of  oscillatory  combustion  instability  in  solid 
rocket  motors  results  from  the  responsiveness  of  the  combustion  pro¬ 
cess  to  oscillations  in  the  flow  environment.  Because  of  the  high 
probability  of  combustion  instability  in  low  signature  motors  and  the 
attendant  jeopardy  to  successful  motor  and  system  performance,  the 
capability  of  instability  prediction  is  of  unquestioned  importance. 

The  state-of-the-art  in  the  field  of  linear  analysis  of  combustion 
instability  is  based  on  a  perturbation  of  the  acoustic  field  in  the 
burning  propellant  and  an  evaluation  of  the  growth/decay  coefficient 
associated  with  the  acoustic  pressure  waves;  a  positive  coefficient 
indicates  an  amplification  of  the  waves  and,  therefore,  instability, 
and  a  negative  value  implies  attenuation  of  the  waves  and  stability. 

In  early  1975  the  development  of  a  Standardized  Stability 
Prediction  Method  for  Solid  Rocket  Motors  [1]  was  begun  by  Aerojet 
Solid  Propulsion  Company  under  contract  with  the  Air  Force.  It  was 
decided  that  this  method  would  employ  the  NASTRAN  finite  element 
fluid  analysis  option  which  was  developed  for  rocket  acoustic  cavity 
analysis  [2,3,4 ,5] ^  The  NASTRAN  finite  element  analysis  is  axisym- 
metric  and  approximates  the  effect  of  radial  slots  (Figure  1)  on  the 
cavity  acoustics.  As  long  as  the  slot  width  is  relatively  narrow,  the 
approximation  provides  an  accurate  model  of  the  cavity  acoustic 
response,  but  as  the  slot  width  increases,  the  accuracy  diminishes. 

The  NASTRAN  option  as  used  in  the  standardized  method  does  not  provide 
for  the  coupling  of  the  vibratory  response  of  the  solid  propellant 
grain  to  that  of  the  acoustic  cavity  and  hence  does  not  provide  a 
means  of  evaluating  the  damping  of  the  acoustic  oscilliations  by  the 
propellant  grain  response.  An  additional  limiting  feature  of  the 
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standardized  method  is  that  it  employs  a  post  processor,  separate 
from  the  other  analyses,  for  evaluating  the  growth/decay  coefficient. 


A< 


Figure  1.  Three-dimensional  cavity-solid  propellant  model. 

The  three-dimensional  analysis  package  presented  herein  was 
developed  primarily  to  provide  1)  more  generality,  and  therefore 
more  accuracy  in  modeling  complex  cavity  geometries,  2)  a  means  of 
predicting  the  damping  of  acoustic  oscillations  by  the  solid  propel¬ 
lant  grain,  3)  an  integrated  program  designed  solely  for  the  purpose 
of  combustion  instability  prediction  with  ease  of  use,  and  4)  a  means 
of  coupling  the  combustion  instability  analysis  to  an  internal  bal¬ 
listics  analysis.  The  developed  package  possesses  all  of  these 
positive  features  as  well  as  certain  inherent,  somewhat  unappealing 
features  associated  with  "bigness"  which  will  be  pointed  out  later 
in  the  paper. 

The  main  aspects  of  the  developed  three-dimensional  code  are 
listed  as  follows  and  will  be  elaborated  upon  subsequently  in  the  paper: 

1)  It  utilizes  a  three-dimensional  finite  element  mesh 
generator  especially  adapted  to  provide  input  to  the 
program. 

2)  It  utilizes  the  principle  of  dihedral  symmetry  which  enables 
a  consideration  of  only  the  smallest  repeating  geometrical 
segment. 
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3)  It  couples  the  response  of  the  gas  cavity  with  that  of  the 
solid  propellant  grain  to  enable  the  calculation  of  the 
frequency  of  the  coupled  system  and  the  damping  provided 
by  the  propellant  grain. 

4)  It  provides  for  modeling  the  propellant  grain  as  a  nearly 
incompressible  material  (which  differs  from  the  common 
minimum  potential  energy  formulation) . 

5)  It  utilizes  the  principle  of  condensation,  wherein,  in 
this  case,  the  fluid  pressure  degrees-of-freedom  are 
designated  ,,master,,  and  the  solid  propellant  displace¬ 
ment  and  mean  pressure  parameter  degrees-of-freedom  are 
designated  "slave."  This  enables  a  major  reduction  in  the 
size  of  the  problem;  the  number  of  equations  is  reduced 
from  the  total  number  of  degrees-of-freedom  of  the  coupled 
system  to  the  number  of  fluid  pressure  degrees-of-free- 
dom. 

6)  It  provides  the  option  of  considering  the  response  of  the 
gas  cavity  alone  (which  models  the  cavity  boundaries  as 
being  "acoustically  hard") .  This  option  might  be  utilized 
in  certain  cases  where  a  savings  in  computer  costs  or 
storage  is  a  dominant  consideration.  In  this  case  the 
condensation  approach  previously  described  obviously 
would  not  be  taken. 

7)  It  calculates  the  three-dimensional  potential  flow  field. 

8)  It  evaluates  the  stability  integrals  for  the  calculation 
of  the  net  driving/damping  coefficient  for  each  mode. 

FINITE  ELEMENT  MESH  GENERATION 

It  is  only  necessary  to  develop  a  finite  element  mesh  for  one 
repeating  segment  (Figure  1)  of  the  total  cavity-grain  rocket  geometry. 
This  is  true  because  of  the  employment  of  the  principle  of  dihedral 
symmetry  in  the  program.  Although  the  three-dimensional  element  used 
in  the  program  for  both  the  cavity  region  and  the  propellant  grain 
is  a  tetrahedron  [6],  the  mesh  is  that  of  bricks,  each  brick,  or 
quasi-hexahedron,  being  comprised  of  five  basic  tetrahedra.  The  break¬ 
down  of  the  quasi-hexahedra  into  tetrahedra  is  performed  internally. 

The  mesh  generator  is  an  efficient  routine  which  automatically  creates 
the  complete  finite  element  mesh,  for  both  cavity  and  solid  pro¬ 
pellant  regions,  from  a  minimal  amount  of  input.  Each  repeating 
segment  is  sectioned  in  the  longitudinal  or  z-direction,  with  each 
section  comprised  of  a  number  of  quadrilateral  parts  which  are 
identified  by  a  counterclockwise  listing  of  their  part  boundary  curves. 
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Part  boundary  curves  may  be  ellipses  as  well  as  straight  lines  and 
their  points  of  intersection  are  designated  by  I,J  indices  as  shown 
in  Figure  2. 


! 


IS,  Jt 


PART  1  (CAVITY  I 

II.  J1  TOI6.J1  (CURVE  1) 

IS,  J1  TO  IS.  Jt  (CURVE  21 
IS,  Jt  TO  11.  Jt  (CURVE  3) 

II.JtTOII.JI  (CURVE  4) 

PART  2  (CAVITY  SOLID  PROPELLANT  INTERFACE) 

II,  J4TOI6.  J4  (CURVE  31 

IS.  Jt  TO  IS,  JS  (CURVE  21 

IS.  JSTOM.  JS  (CURVE  SI 

II.  JS  TO  II.  Jt  (CURVE  41 
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Figure  2.  Finite  element  mesh  definition  for  cavity-solid 
propellant. 

By  coupling  the  output  of  a  standard  internal  ballistics 
analysis  computer  code  with  the  mesh  generator  one  can  obtain  a 
combustion  instability  analysis  at  any  stage  of  performance,  i.e.,  at 
any  instant  from  ignition  of  the  propellant  to  burnout.  This  has 
been  done  through  the  definition  of  additional  part  boundary  curves, 
which  conform  to  the  geometry  of  the  different  zones  of  burning  for 
the  differing  cavity-grain  configurations.  With  the  transfer  of  a 
small  number  of  geometric  parameters  from  the  output  of  the  internal 
ballistics  code  to  the  mesh  generator,  a  complete  combustion  insta¬ 
bility  analysis  can  be  initiated  and  automatically  performed.  The 
designation  of  these  parameters  and  a  detailed  description  of  the 
development  of  the  coupling  program  are  found  in  References  [7]  and  [8]. 
Presently  the  coupling  of  an  internal  ballistics  analysis  with  the 
combustion  instability  analysis  is  limited  to  the  consideration  of 
axisymmetric  and  star  cross-section  geometries,  but  increased  part 
boundary  curve  definition  is  a  matter  of  expansion  of  the  program, 
not  of  additional  development. 
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DIHEDRAL  SYMMETRY 

If  a  geometrically  defined  body  is  made  up  of  identical  seg¬ 
ments  symmetrically  arranged  with  respect  to  an  axis,  the  degrees- 
of-freedom  for  a  finite  element  analysis  can  be  transformed  into 
uncoupled  symmetrical  components,  thereby  greatly  reducing  the  number 
of  equations  which  must  be  solved  simultaneously  [9].  A  further 
reduction  can  be  effected  if  each  segment  has  a  plane  of  reflective 
symmetry.  Dihedral  symmetry  is  the  term  applied  to  this  latter 
condition.  It  can  be  seen  from  Figure  1  that  a  typical  solid  rocket 
geometry  meets  this  requirement  and,  therefore,  the  principle  of 
dihedral  symmetry  can  be  employed  in  a  three-dimensional  (cylindrical 
coordinate)  analysis.  The  application  of  the  principle  to  this  prob¬ 
lem  is  explained  in  detail  in  Reference  [10]  and  will  not  be  reproduced 
here,  but  the  resulting  analysis  will  be  discussed. 

Employing  the  principle  of  dihedral  symmetry,  the  program 
computes  three  distinct  types  of  acoustic  harmonics:  the  zero 
harmonic,  the  K  harmonics,  and  the  N/2  harmonic.  The  zero  harmonic 
exists  for  all  cases.  The  number  of  possible  K  harmonics  is  given  by 

K  =  1,  .  .  . ,  L  (1) 

where  the  harmonic  index  L  is  given  by 


(if  N  is  an  odd  number) 


(2a) 


(if  N  is  an  even  number) 


(2b) 


where  N  is  the  number  of  radial  slots  or  lobes.  The  N/2  harmonic 
exists  only  when  N  is  even.  Referring  to  the  geometry  of  Figure  1, 
then,  one  could  calculate  the  zero,  first,  and  second  harmonics,  the 
latter  two  being  K  harmonics.  The  longitudinal  modes  associated  with 
each  harmonic  are  calculated  as  requested  by  the  user.  The  result 
of  this  calculation  (eigensolution)  is  the  natural  circular  frequency 
associated  with  each  acoustic  mode  and  the  corresponding  acoustic 
pressure  distribution  (normalized  acoustic  pressure  at  each  finite 
element  nodal  point)  for  the  smallest  repeating  segment.  The  pressure 
distribution  for  each  of  the  other  segments  is  then  simply  calculated 
automatically  through  the  dihedral  transformation  matrix.  The 
acoustic  velocity  components  (constant  for  the  region  occupied  by  each 
cavity  tetrahedral  element)  are  computed  from  the  acoustic  pressure  nodal 
point  values  for  those  cavity  elements  which  are  adjacent  to  the  propellant 
grain. 
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The  theoretical  formulation  of  the  complete  three-dimensional 
finite  element  analysis,  in  which  the  natural  circular  frequency,  the 
acoustic  pressure  distribution,  and  the  element  acoustic  velocities 
are  calculated,  is  given  in  Reference  [11]  and  will  not  be  repeated 
here.  Reference  [12]  also  presents  the  theoretical  finite  element 
formulation  which  was  used  in  the  development  of  this  program. 


COUPLED  RESPONSE 

The  presence  of  the  solid  propellant  grain  can  significantly 
shift  the  acoustic  system  frequency  from  that  of  the  gas  phase  alone, 
a  portion  of  the  acoustic  energy  being  dissipated  by  the  deformable 
solid  material.  This  effect  can  be  one  of  the  more  significant  sinks 
for  acoustic  energy  in  both  large  and  small  rocket  motors,  the  amount 
of  damping  depending  on  the  grain  geometry  and  mechanical  properties 
and  on  the  acoustic  mode  shape  and  natural  frequency. 

In  order  to  evaluate  the  coupled  cavity-solid  propellant 
grain  response  it  is  necessary  to  model,  by  the  finite  element  method, 
both  cavity  and  grain.  This,  of  course,  greatly  increases  the  size  of 
the  problem  to  be  solved,  from  the  standpoint  of  number  of  initial 
degrees-of-freedom.  The  coupled  finite  element  formulation  of  Refer¬ 
ence  [12],  which  is  coded  in  the  three-dimensional  program,  is 
expressed  in  matrix  form  by: 
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where  [F]  is  the  fluid  inertia  matrix,  [C ]  is  the  fluid  compressibil¬ 
ity  matrix,  [K]  is  the  solid  stiffness  matrix,  [M]  is  the  solid  con¬ 
sistent  mass  matrix,  [U]  is  the  matrix  which  couples  acoustic  pressure 
degrees-of-freedom  to  solid  displacement  degrees-of-freedom,  { p }  is 
the  acoustic  pressure  vector,  {A}  is  the  solid  displacement  vector, 

{H}  is  the  mean  pressure  parameter  vector  (to  be  elaborated  upon  later), 
and  r  is  the  eigenvalue  of  the  coupled  system. 

The  structural  damping  can  be  attributed  to  the  out-of-phase 
response  of  the  solid  propellant  grain,  which  is  measured  in  terms  of 
the  complex  shearing  modulus  of  the  grain,  which,  in  turn,  results  in 
a  complex  eigenvalue  for  the  coupled  system.  The  imaginary  part  of 
the  complex  eigenvalue  obtained  from  the  eigensolution  is  the  natural 
circular  frequency  of  the  coupled  system  while  the  real  part  is  the 
structural  damping  rate. 


HACKETT 


Since  the  complex  shearing  modulus  is  frequency  dependent,  a 
series  of  iterations  may  be  necessary  before  the  accurate  value  of 
complex  modulus  for  input  into  the  program  is  determined. 


PROPELLANT  GRAIN  MODELING 

Since  the  propellant  grain  is  only  accurately  modeled  as  a 
nearly  incompressible  material,  the  well-known  standard  Navier  dis¬ 
placement  formulation,  in  conjunction  with  the  Ritz  procedure,  would 
lead  to  inaccuracies  in  the  finite  element  modeling  of  the  grain.  In 
order  to  avoid  this  situation  the  solid  finite  element  formulation 
utilized  in  this  program  is  that  of  a  linear  displacement-linear  mean 
pressure  tetrahedron  [13].  It  is  similar  to  the  Herrmann  variational 
formulation  [14],  which  employs  a  linear  displacement  function  and  a 
constant  mean  pressure  function.  The  finite  element  modeling  of  the 
propellant  grain  used  in  this  program  is  outlined  in  detail  in 
Reference  [12]  and  will  not  be  repeated  here. 


EIGENVALUE  ECONOMIZER  -  CONDENSATION 


The  extraction  of  eigenvalues  and  eigenvectors  is  a  much  more 
expensive  operation  than  is  the  solution  of  simultaneous  linear 
equations.  It  requires  roughly  twice  as  nueh  time  to  extract  a  single 
eigenvalue  as  to  do  a  single  "static"  analysis.  In  order  to  reduce  or 
"condense"  the  number  of  degrees-of-f reedom  in  the  eigensolution ,  the 
following  technique  is  utilized  in  the  three-dimensional  program.  Fur¬ 
ther  details  of  the  method  are  found  in  Reference  [15]. 

The  original  formulation  of  the  coupled  system  is  given  by 
Equation  3,  where  the  number  of  degrees-of-freedom  is  equal  to  the  num¬ 
ber  of  cavity  nodal  point  pressures,  plus  the  number  of  solid  propellant 
nodal  point  displacement  components,  plus  the  number  of  solid  propellant 
nodal  point  mean  pressure  parameter  values  (one  at  each  propellant  grain 
nodal  point),  for  the  analyzed  repeating  segment.  The  condensed  formu¬ 
lation  is  given  by: 


<[Fr]  '  '2  [Cr]  >  |  -)  ’  0 
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and 


(5b) 


where  I  is  the  identity  matrix.  The  relationship  between  initial  and 
reduced  degrees-of- freedom  is 


UNCOUPLED  RESPONSE 


Although  one  of  the  most  important  features  of  the  three- 
dimensional  program  is  the  coupled  cavity-propellant  response,  the 
option  of  a  cavity  analysis  alone  is  available  to  the  user.  In  this 
case  the  cavity-solid  propellant  interface  is  modeled  as  an  "acousti¬ 
cally  hard"  boundary.  The  cavity  only  option  would  greatly  reduce 
the  solution  time  and  cost  and  would,  in  some  cases,  perhaps  suffice. 
The  user  need  make  only  two  simple  adjustments  to  the  input  data,  and 
these  adjustments  are  described  in  the  user's  guide  to  the  three- 
dimensional  code  [16]. 


POTENTIAL  FLOW  CALCULATION 


A  separate  operation  is  carried  out  in  the  sub-routine  which 
performs  a  potential  flow  analysis  for  the  purpose  of  determining  the 
mean  flow  field  in  the  rocket  cavity*  As  in  the  case  of  the  eigen- 
solution,  only  the  smallest  repeating  geometrical  segment  need  be 
considered,  and,  for  this  calculation,  only  the  cavity  portion  of  the 
segment  with  the  proper  boundary  conditions.  The  same  general  formula¬ 
tion  of  the  finite  element  model  equations  of  motion  is  utilized  except 
that,  in  this  case,  the  fluid  is  considered  to  be  incompressible.  The 
mass  flow  into  the  cavity  from  the  burning  propellant  surface  is 
modeled  as  a  cavity-solid  propellant  interface  nodal  point  quantity. 

It  is  calculated  by  summing  the  interface  surface  areas  associated 
with  each  nodal  point  lying  on  the  cavity-solid  propellant  interface. 
The  solution  of  the  resulting  set  of  linear  equations  for  the  mean 
flow  velocity  components  (constant  for  each  cavity  tetrahedron)  is 
explained  in  Reference  [11]  and  will  not  be  discussed  here. 
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EVALUATION  OF  STABILITY  INTEGRALS 

The  final  operation  consists  primarily  of  the  evaluation  of 
the  stability  integrals  associated  with  the  various  driving/ damping 
coupling  mechanisms  which  occur  in  the  cavity  chamber  in  the  presence 
of  combustion  and  flow*  The  stability  integrals  presently  incorporated 
into  the  program  are  those  derived  in  Reference  [17]  for  the  three- 
dimensional  case,  along  with  the  flow-turning  formulation  also  found 
in  Reference  [17]*  The  use  of  linear  pressure  (and  therefore  constant 
velocity)  tetrahedral  elements  to  represent  the  cavity  region  enables 
an  explicit  evaluation  of  the  stability  integrals,  given  the  acoustic 
nodal  point  pressures  and  element  velocities  and  the  mean  flow  element 
velocities  from  the  finite  element  solutions.  At  present  the  code 
does  not  contain  a  routine  for  evaluating  nozzle  damping  or  particu¬ 
late  damping.  These  additional  calculations  will  be  added. 

The  evaluation  of  stability  integrals  is  for  the  purpose  of 
calculating  the  driving/damping  coefficient,  a,  a  fact  well-known  to 
the  combustion  community.  A  positive  a  indicates  a  stable  mode.  The 
net  value  of  a  computed  by  the  three-dimensional  code  is  a  summation 
of  the  computed  values  of  apc  (pressure  coupling),  a^c  (velocity 

coupling),  aFT  (flow-turning),  and  agD  (structural  damping).  The  value 

of  is  obtained  from  the  complex  eigensolution  described  in  an 

earlier  section;  the  other  three  a-values  are  obtained  from  the  evalu¬ 
ation  of  the  stability  integrals.  It  is  known  that  the  pressure 
coupling  mechanism  always  drives  the  acoustic  oscillations,  that  the 
velocity  coupling  mechanism  may  either  drive  or  damp  the  oscillations, 
and  that  the  flow-turning  mechanism  always  damps  the  oscillations. 
Response  factors  [18,19]  are  input  into  the  program  as  multiples  of  the 
stability  integrals  for  the  calculation  of  the  afs  obtained  from  the  dif¬ 
ferent  coupling  mechanisms.  These  propellant  grain-dependent  response 
functions  are  obtained  from  other  analyses  and  utilized  as  direct 
input  into  the  three-dimensional  program. 

COMPARISON  WITH  ANALYTICAL  SOLUTIONS  AND 
EXPERIMENTAL  OBSERVATIONS 

In  order  to  affirm  the  accuracy  and  usefulness  of  any  numer¬ 
ical  analysis  package  it  is  necessary  to  make  comparisons  of  results 
obtained  numerically  with  available  analytical  closed-form  solutions 
and  experimentally  obtained  results.  In  the  case  of  a  complex  anal¬ 
ysis  procedure  such  as  that  of  combustion  instability,  available  bases 
of  comparison  are  limited  in  both  regards.  The  following  results  are 
presented  to  support  the  contention  of  accuracy  and  utility  of  the 
developed  package. 
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Analytical  Comparisons.  A  right  circular  cylindrical  cavity 
254mm  long  and  254mm  in  diameter  was  modeled  using  an  18°  repeating 
segment  consisting  of  32  quasi-hexahedral  elements,  4  in  the  radial 
direction  and  8  in  the  longitudinal  direction.  Acoustic  frequencies 
corresponding  to  the  first  and  second  logitudinal  modes  (1L  and  2L) , 
the  first  tangential  mode  (IT),  the  first  tangential-first  longitudinal 
mode  (1T-1L) ,  and  the  first  radial  mode  (1R)  were  solved  for  and  are 
compared  with  the  corresponding  analytically  obtained  frequencies  of  the 
cavity.  The  values  are  shown  in  Table  1. 

TABLE  1.  NUMERICALLY  AND  ANALYTICALLY  CALCULATED  ACOUSTIC 
FREQUENCIES  FOR  CYLINDRICAL  CAVITY 


Frequency  (Hz) 

Mode 

Analytical 

Solution 

Finite  Element 
Solution 

1L 

2000 

1988 

2L 

4000 

4042 

IT 

2400 

2378 

1T-1L 

3124 

3151 

1R 

4994 

5033 

As  described  earlier,  a  separate  potential  flow  calculation 
is  made  in  the  combustion  instability  analysis  package  for  the  purpose 
of  evaluating  the  stability  integrals.  A  comparison  of  the  finite 
element  potential  flow  solution  with  the  analytical  solution  for  the 
cylindrical  cavity  is  shown  in  Figure  3. 


y 


Figure  3.  Comparison  of  finite  element  and  analytical  potential 
flow  field  solutions  for  a  circular  cylinder. 
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The  accuracy  of  the  finite  element  frequency  calculations  is 
seen  in  Table  1.  These  results  are  highly  favorable,  especially  when 
considering  the  relatively  small  number  of  elements  used  to  model  the 
cavity.  The  somewhat  less  satisfying  agreement  in  the  case  of  the 
finite  element  potential  flow  solution  can  be  attributed  to  the 
relatively  small  number  of  elements,  especially  near  the  exhaust  end 
of  the  cavity  where  the  axial  velocity  component  has  its  largest  value. 
Additional  comparisons  with  closed-form  solutions  of  the  same  type  have 
yielded  results  similar  to  those  cited  here. 

Experimental  Comparisons.  The  experimental  data  used  for 
comparison  purposes  were  obtained  from  Hercules,  Inc.,  Allegany 
Ballistics  Laboratory  (ABL)  in  Cumberland  Maryland  and  result  from 
testing  of  an  early  experimental  model  of  the  Smokeless  CHAPARRAL 
Motor  Assembly,  conducted  in  early  1977.  The  finocvl  geometry  of 
the  pulsed  motor  and  the  finite  element  model  consisted  of  a  small 
diame ter-to-length  ratio  cylindrical  cavity  having  four  radial  slots 
over  approximately  one-third  of  the  cavity  length.  Experimental 
longitudinal  mode  frequency  and  stability  data  can  be  compared  to 
corresponding  finite  element  analysis  results  in  Table  2.  The  test 
data  and  numerical  results  found  in  Table  2  correspond  to  a  burn 
distance  of  8.03mm.  The  finite  element  analyses  were  performed  on 
the  CDC  6600  computer  at  Redstone  Arsenal.  The  lack  of  close  agree¬ 
ment  observed  when  comparing  corresponding  net  a-values  can  to  some 
extent  be  attributed  to  the  use  of  somewhat  inaccurate  values  of 
response  factors,  which  are  presently  considered  in  the  program  as 
properties  of  the  propellant,  but  which  are  also  flow-dependent. 
However,  the  qualitatively  good  agreement  between  corresponding  values 
can  be  observed. 

TABLE  2.  SMOKELESS  CHAPARRAL  DATA  AND  FINITE  ELEMENT  ANALYSIS  RESULTS 


Test 

Data 

Finite  Element  Analysis 

Frequency  (Hz) 

a  (sec~^) 

Frequency  (Hz) 

a  (sec~l) 

323 

— 

-94 

290 

-41 

307 

-198 

644 

-266 

640 

-290 

593 

-220 

909 

-192 

927 

t 

-145 

909 

-148 
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The  results  from  additional  finite  element  analyses  of  the 
early  experimental  model  of  the  Smokeless  CHAPARRAL  Motor  Assembly  are 
shown  in  Table  3.  The  most  notable  aspect  of  this  set  of  finite  element 
analysis  results  is  the  prediction  of  instability  in  some  of  the 
tangential  modes  (positive  values  of  a) .  This  compares  quite  favorably 
with  ABL  tests  of  the  same  assembly,  which  exhibited  clear  instabilities 
in  tangential  modes.  Also,  ABL  test  results  showed  5-6  separate  tangen¬ 
tial  modes  evident  in  the  frequency  range  of  9000-20,000  Hz. 

TABLE  3.  HIGHER  ACOUSTIC  MODE  FINITE  ELEMENT  ANALYSIS  RESULTS 
FOR  SMOKELESS  CHAPARRAL  MOTOR 


Harmonic 

Mode 

Frequency 

(Hz) 

a  (sec~^) 

1 

7992 

-147 

2 

13,476 

+  33 

K  (=1) 

3 

15,303 

-  55 

4 

15,915 

+  46 

_ 

5 

17,907 

-  92 

1 

8981 

-232 

! 

2  1 

15,948 

+162 

N/2  (=2) 

3 

20,516 

-119 

4 

23,693 

-  62 

5 

25,644 

-  85 

The  comparisons  between  numerical  and  experimental  results, 
which  can  be  made  from  Tables  2  and  3,  would  appear  to  support  the 
contention  that  the  developed  finite  element  package  can  be  of 
valuable  use  in  predicting  the  instability  of  solid  rocket  systems. 

CONCLUSIONS 

The  developed  computer  package  performs  a  linear  analysis  of 
the  irrotational  motions  of  an  inviscid,  compressible  fluid  coupled 
to  the  motions  of  a  nearly  incompressible,  linearly  viscoelastic  solid, 
performs  a  linear  potential  flow  analysis  of  the  irrotational  motions 
of  an  inviscid,  incompressible  fluid,  and  then  determines  the  effect 
of  the  flow  field  and  of  combustion  on  the  calculated  acoustic 
oscillations.  There  are  obvious  limitations  attached  to  any  code 
which  is  as  basic  as  the  restrictions  listed  above  dictate,  but  it  is 
felt  that  the  developed  code  presented  herein  is  probably  as  sophis¬ 
ticated  as  the  present  state-of-the-art  warrants.  It  is  viewed  as 
having  much  potential  as  both  a  design  and  a  research  tool.  As  the 
state-of-the-art  in  combustion  technology  advances,  it  is  felt  that 
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the  code  can  be  relatively  easily  revised  and  updated  to  include  the 
new  technology;  at  least  it  was  designed  with  that  provision  in 
mind . 

One  of  the  most  attractive  features  of  the  code  is  the  ease 
of  use.  Other  extremely  important  attributes  are  the  fact  that  it 
is  three-dimensional,  that  it  performs  a  coupled  cavity-solid  propel¬ 
lant  analysis  (or,  alternatively,  a  cavity-only  analysis),  and  that  all 
analyses  are  automatically  linked.  Features  of  the  code  which  do  not 
enhance  its  reputation  also  exist,  and  they  too  should  be  pointed  out. 
It  is  a  large  program  requiring  a  large  amount  of  storage  and  it  may 
require  long  run  times,  as  is  the  case  with  any  three-dimensional 
finite  element  program.  Presently,  the  entire  program  is  in-core 
computation,  but  this  will,  in  all  probability,  be  modified.  The 
eigenvalue  routine  used  in  the  analysis  may  not  be  the  most  efficient 
one  available.  This  possibility  is  currently  being  investigated,  and 
if  a  more  efficient  routine  can  be  found,  it  will  be  used  in  place  of 
the  one  presently  employed.  In  certain  instances  a  two-dimensional 
(axisymmetr ic)  uncoupled  analysis  provides  sufficient  accuracy,  and 
for  such  cases  use  of  the  three-dimensional  code  might  not  have  merit. 

It  is  felt  that  the  demonstrated  attributes  of  the  three- 
dimensional  code  far  outweigh  any  foreseen  disadvantages,  and  that 
it  can  provide  the  means  of  performing  important,  and  heretofore 
impossible,  analyses. 
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